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Abstract 



We investigate the experimental implications of the minimal gauge- 
mediated low energy supersymmetry breaking (GMLESB) model for Fermilab 
Tevatron collider experiments. We map out the regions of parameter space 
of this model that have already been excluded by collider searches and by 
limits on 6 ^ 57. We use ISAJET to compute the cross sections for a variety 
of topological signatures which include photons in assocation with multiple 
leptons, jets and missing transverse energy. The reach in the parameter A, 
which fixes the scale of sparticle masses, is estimated to be ~ 60, 100 and 135 
TeV for Tevatron integrated luminosities of 0.1, 2 and 25 fb~^, respectively. 
The largest signals occur in photon(s) plus lepton(s) plus multi-jet channels; 
jet-free channels containing just photons plus leptons occur at much smaller 
rates, at least within this minimal framework. 
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I. INTRODUCTION 



The search for weak scale supersymmetry (SUSY) |Q forms an integral part of the ex- 
perimental program at all high energy colliders in operation or in the construction and 
planning phases. In the absence of any one compelling theoretical framework, the experi- 
mental analyses have to be performed within the context of particular models. The Minimal 
Supersymmetric Standard Model (MSSM) is obtained by the direct supersymmetrization of 
the Standard Model (SM), but with two Higgs doublets, and including all renormalizable soft 
supersymmetry breaking interactions consistent with SM symmetries. The resulting theory 
has over one hundred model parameters making phenomenological analyses intractable. (If 
i?-parity violation by renormalizable baryon- (or lepton-) number violating superpotential 
interactions is allowed, the number of model parameters is even larger.) The proliferation of 
soft-SUSY breaking parameters is a reflection of our lack of understanding of the mechanism 
of SUSY breaking. The practical solution for reducing the parameter-space is to incorporate 
simplifying ansatze usually based on the assumed symmetries of physics at very high scales. 

One especially attractive and economic realization of this idea is provided by the so- 
called minimal supergravity (mSUGRA) framework that has been recently used for most 
phenomenological 0, and also some experimental IQ, analyses of SUSY. Here, "minimal" 
refers in part to the technical assumption of canonical kinetic energy terms. It is envisioned 
that SUSY is dynamically broken at a scale Msusy in a sector of the theory that interacts 
with the observable sector of quarks, leptons, gauge and Higgs bosons and their superpart- 
ners only via gravity, which acts as the "messenger" of supersymmetry breaking As a 
result the particle-sparticle mass gap in the observable sector is suppressed by relative 
to Msusy and is thus given by ~ Mg^gy/Mp. This quantity may be of order the weak 
scale if Msusy is ~ 10^^ GeV, where the reduced Planck mass Mp = 2.4 x 10^^ GeV. The 
Goldstone fermion (which dominantly lives in the hidden sector) then forms the longitudinal 
components of the gravitino which generically acquires a mass ~ M^jj^y/^P by the super- 
Higgs mechanism 0. Although a weak scale particle, the couplings of the gravitino are of 
gravitational strength, so that it plays no role in particle physics. The resulting low energy 
Lagrangian ^ in the observable sector is just a globally supersymmetric Lagrangian with 
universal scalar (mo) and gaugino masses {mi/2) and a universal trilinear scalar soft-SUSY 
breaking parameter {Aq) at an ultra-high scale Mx often identified with Mqut- The univer- 
sality of the gauginos may have its origins in grand unification while the universal boundary 
condition for the scalar masses results from our technical assumption of the canonical ki- 
netic energy terms, mentioned above. Although these boundary conditions are not generic 
1^ to supergravity models (and are tantamount to assuming an additional global symmetry 
known to be broken by Yukawa interactions) this framework is generally referred to as the 
minimal SUGRA framework. A very attractive feature of this picture is that over a signifi- 
cant portion of the parameter space of the model, radiative corrections lead to the correct 
pattern of electroweak symmetry breaking the SUSY Higgs mass parameter /x^ is then 
fixed by the value of Mz- In such a scenario, all the sparticle properties are determined 
by just four additional parameters along with sgn /i. The lightest supersymmetric particle 
(LSP) is frequently the lightest neutralino {Zi) and is a good candidate for cosmological 
cold dark matter if i?-parity is conserved as is assumed to be the case 0. 

Attractive though this framework is, it says nothing about the dynamics of SUSY break- 
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ing. In recent studies, Dine, Nelson, Nir and Shirman have attempted to construct 
models where SUSY is dynamically broken in the hidden sector of the theory and commu- 
nicated to a messenger sector via new gauge interactions. The messenger sector, which is 
characterized by a mass scale M, interacts with the observable sector via the known SM 
gauge interactions which then serve to communicate SUSY breaking to the visible sector 
quarks, leptons, gauge and Higgs bosons, and their superpartners. The effective SUSY 
breaking scale in the observable sector is now suppressed by M rather than Mp and is 
~ ^ X Mgjjgy/^y with MsusY being the induced SUSY breaking scale in the messenger 
sector and a is the relevant SM fine structure constant. The effective scale of SUSY break- 
ing in the observable sector may thus be ~ Mw^^k even if the SUSY breaking scale and the 
messenger scale M are as small as few tens or few hundred TeV. The gravitino mass, which 
is still suppressed by Mp, is smaller by a factor M/Mp ~ 10~^^ compared to rriQ within the 
SUGRA framework for M ~ 250 TeV. Thus in these new scenarios, the gravitino mass may 



be in the electron-volt range. Fayet [11 1 has shown that for the longitudinal components 



of such a superlight gravitino the smallness of the gravitational coupling is made up by the 
size of the wavefunction of the gravitino of electroweak scale energy so that this gravitino 
does not decouple from other particles. 

The fact that the gravitino, and not the lightest neutralino, is the LSP is the main reason 



why the phenomenology of these models can be quite different |12|-[T^ from the usual MSSM 
analyses. Most importantly, the lightest neutralino (as well as other neutralinos) can now 
decay via Zi 6*7, and also via Zi GZ or GHi (where i = i, h or p for light, heavy 
and pseudoscalar Higgs bosons) if the decays are kinematically allowed. We will see that 
for sparticle masses accessible at the Tevatron, the photon decay dominates over much of 
the parameter space of the model. Since the gravitino escapes experimental detection, we 
expect that in such a scenario SUSY events will generically have the n-jet(s)-|-m-lepton(s)-|-fc- 
7 + topology. We find that the branching ratio for sparticles other than Zi to decay via 
the gravitino mode is small, so that k = — 2 because the photon detection efficiency is not 
unity. This novel source of photons in SUSY events was considered |[T^-[1^ to be the origin 
of the single e"'"e~77 -|- event |15| recorded by the CDF Collaboration. Finally, we will 
show that at least within the simplest of the gauge mediated low energy SUSY breaking 
(GMLESB) models reviewed in the next Section, there would have been a plethora of other 
events accompanying the CDF event, making its SUSY origin within this context rather 
implausible. 

The scenarios envisioned in Ref. |TD[ are in a sense considerably more ambitious than the 
conventional SUGRA picture since they include not only a mechanism for the transmission of 
SUSY breaking, but also the dynamical mechanism for it. Thus, all scales in the low energy 
theory, in particular the values of fi and the SUSY breaking Higgs boson mass parameters 
that describes the observable sector, should be derived in such a scenario. Since there is 
no universally accepted resolution of the fi problem |T6[, we will adopt a phenomenological 
approach and focus on the implications of the GMLESB model, treating /i to be a parameter 
that gets fixed by the value of Mz via the constraint from radiative electroweak symmetry 



breaking |1J] . In other words, we regard the mediation of SUSY breaking and the mechanism 
of SUSY breaking as independent issues, with independent consequences. 

The remainder of this paper is organized as follows. In Section II we briefly review 
the assumptions underlying the GMLESB model and set up the parameter space for our 



3 



phenomenological analysis. We delineate the regions of parameter space already excluded 
by experiments at LEP2 and the measurement of the 6 — >■ 57 branching ratio by the CLEO 
Collaboration. Finally, in Sec. II we study the branching fractions for the direct decays of 
sparticles to gravitinos which is the novel feature of these scenarios. In Section III, we use 
ISAJET to generate events which lead to the various n— j + m — i + k — '-/ + I^j, event 
topologies at the Fermilab Tevatron and give an estimate of its reach in various channels 
from the data of Run I as well as from Run II with the Main Injector and the proposed 
TeV33 upgrade. We summarize our results in Section IV. 



II. THE MINIMAL MODEL OF GAUGE MEDIATED SUPERSYMMETRY 

BREAKING 

A. Model Parameter Space 

Supersymmetry is dynamically broken in the "secluded sector" (this was referred to as the 
hidden sector in the SUGRA framework) of the theory and communicated to the "messenger 
sector" via some new gauge interactions which do not couple to the known particles. In the 
simplest realization of this idea |TD[ the messenger sector is weakly coupled (so that non- 



perturbative dynamics does not cause SUSY breaking via gaugino masses) and comprises of 
one set of "quark" and "lepton" superfields in a 5 + 5 representation of SU{^) coupled to 

a singlet via a superpotential of the form W = XiSqq + \2Sii. The incorporation of new 
fields in complete GUT multiplets ensures that they do not spoil the successful prediction of 
sin^ 6\Y if this model is incorporated into a GUT. The scalar and auxilliary components of the 
field S acquire vacuum expectation values, {S) and (F), the latter signalling the breaking of 
SUSY in the messenger sector. SM gauge interactions then carry the information of SUSY 
breaking to the observable sector, and induce masses (proportional to the corresponding fine 
structure constant) for the gauginos via one loop quantum corrections. The chiral scalars 
feel the effect of SUSY breaking only via these gaugino masses, so that SUSY breaking scalar 
squared masses are induced only as two loop effects. If {F) <^ {S^), the gaugino and scalar 
masses are respectively given by []TU|, 



my = ^A, (2.1) 



and 

,2 _ oa2 



^scalar ~ 2A 



(2.2) 



/2 



with A = {F)/ {S), and ai given in terms of the usual hypercharge coupling g' by ai = . 
Finally, C3 = | for colour triplets and zero for colour singlets while C2 = \ for weak 
doublets and zero for weak singlets. These relations, which are independent of the messenger 
sector superpotential couplings Ai,2, get corrections of ~ (x^^y) which are ignored in the 
subsequent analysis. Notice that instead of a universal scalar mass as in SUGRA, the 
masses of the scalars in this model depend on their gauge quantum numbers: squarks 
are the heaviest, followed by uncoloured electroweak doublets, followed by the colour and 
electroweak singlets. Since the gaugino masses are radiatively generated ||20|, the mass 
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relation is exactly as in a GUT model, although the physics behind this is very different. 
SUSY brealcing A-parameters and the S-parameter are induced only at higher loops so 
that it is reasonable to suppose that these are small. The supersymmetric n parameter is 
not determined by how SUSY breaking is mediated but will be fixed (up to a sign) by the 
constraints of radiative symmetry breaking as in the SUGRA framework. A complete theory 
that includes the dynamics of SUSY breaking will presumably yield a value of /i consistent 
with this. 

Eq. (2.1) and Eq. (2.2) should be regarded as boundary conditions for the gaugino and 
scalar masses valid at the messenger scale M = X{S) (we assume Ai ~ A2), so that these 
parameters need to be evolved down to the weak scale relevant for phenomenological anal- 
ysis. The renormalization group evolution (RGE) of various SUSY breaking gaugino and 
scalar masses is illustrated in Fig. 1, assuming the boundary conditions discussed above. 
In this example we have chosen A = 40 TeV, M = 500 GeV, tan/? = 2 and taken /i < 0. 
The top mass is fixed to mj = 175 GeV throughout this paper. As expected the squarks, 
on account of their QCD interactions, are significantly heavier than all the other scalars. 
The SUSY breaking t-squark masses are smaller than those for other squarks on account 
of their large Yukawa interactions which reduce their masses as Q is evolved down to the 
weak scale. The right-handed sleptons which have only hypercharge gauge interactions are 
considerably lighter than the left-handed slepton or Higgs doublets. Finally, the running 
gaugino masses are proportional to the corresponding fine structure constants at all scales. 
The most important feature of Fig. 1 is that fn'jj^ becomes negative and electroweak sym- 
metry is radiatively broken just as in SUGRA models. Then, we can eliminate the weak 
scale S-parameter in favour of tan (3 (we will return to the issue of whether the resulting 
value of B evolved to the scale M is compatible with the expectation from the boundary 
condition) while //^ is determined by M|. In doing so, we have minimized the one- loop 
corrected effective potential. The model is thus completely specified by the parameter set 
(A, tan/3, M, sgn/x). The dependence on M is presumably logarithmic since it only enters 
via the boundary conditions, so that the A vs. tan (3 plane provides a convenient arena for 
presenting our results. 



B. Sparticle Masses and Experimental Constraints 

We begin by showing contours of various sparticle masses and the weak scale SUSY 
parameters At and fi in the A vs. tan (3 plane for the two signs of /i. We show contours of 
m^^, m^^, and m^^ in Fig. 2a for /i < and Fig. 2b for > 0. In Fig. 2c and Fig. 2d 
we show contours for rriq and rrig in addition to contours for At and /i for the two signs of 
/i. We have fixed M=500 TeV. 

The region in Fig. 2 denoted by bricks is where the proper breaking of electroweak 
symmetry is not obtained. The hatched region is where the lightest neutral Higgs boson 
rriH^ < 60 GeV or < 79 GeV. The latter bound has recently been obtained by the 
ALEPH collaboration at LEP2 [^. The chargino bound is derived assuming that the 
chargino is gaugino-like with ~''^Zi — GeV, and further that the Zi escapes detection. 
Although no analyses have been specifically carried out for the GMLESB scenario, the 
clean experimental environment makes it difficult to imagine that these chargino signals 
would have evaded detection even if Zi were unstable and decayed within the detector via 
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Zx G'y. The LEP limit of ttih > 63 GeV has been obtained for a SM Higgs boson; our 
corresponding requirement of rriH^ > 60 GeV for the hght MSSM Higgs boson should be 
an excellent approximation for the GMLESB framework since the additional Higgs bosons 
are all comparatively heavy (see below). The cross-hatched region at large tan/3 is where 
< m^^ (recall that tl—tr mixing can be substantial for large values of tan /?). While this 
last region would have been excluded by cosmological considerations within the mSUGRA 
framework, this is not so for the GMLESB model since the fi is unstable. Indeed, if the 
model parameters are in the cross-hatched region the phenomenology will be quite different. 
Charginos and neutralinos (including Zi) will cascade decay to fi which will then decay via 
fi — > tG with a width (independent of stau mixing) given by 

m^rG) = -^ ^<-f{ (2.3) 

which yields a decay length of ~ 1.8 x 10^^7fj/5f^(mfj/100G'eK)~^(m(^/lel^)^ cm. Thus, for 
this region of parameters, every SUSY event will contain 2-4 r's in the final state instead of 
hard isolated photons. We will not elaborate this scenario any further. 
The following features of Fig. 2 are worth noting. 

• The scale of sparticle masses is set by A. Except for the cross hatched region, Z\ is the 
next-to-lightest sparticle (NLSP), and further, except for tan/3 close to unity, sparticle 
masses are insensitive to tan /?. 

• The left-handed sleptons (and sneutrinos) are substantially heavier than Ir\ squarks 
are generally heavier than gluinos except for third generation squarks which have 
substantial Yukawa interactions (see Fig. 1). It is worth mentioning that the ratio 

(at the scale M) decreases as the square root of the number of SU (5) vector multiplets 



in the messenger sector [T^ , and can go below unity. If the number of messenger sector 
fields is large, this will effect the phenomenology which is sensitive to whether squarks 
are decaying into gluinos or vice-versa. A similar comment applies to slepton and 
electroweak gaugino masses. 

The value of /i which is obtained from radiative symmetry breaking is large so that 
the lighter chargino and the two lightest neutralinos are gaugino-like, while W2 and 
^3^4 contain large Higgsino components. Thus the chargino and neutralino mass and 
mixing patterns are qualitatively similar to those in the mSUGRA model. 

Even though we start with AtiM^ = 0, an A-parameter of several hundred GeV is gen- 
erated by the renormalization group running: since /i. At and m^^ ^ have comparable 
magnitudes, there is considerable mixing in the t-squark sector. 

Although we have not shown this, all but the lightest of the neutral Higgs scalars 
tend to be heavy. We have checked that these are always heavier than 180 GeV, and 
frequently much heavier. Again, this feature is common with the mSUGRA model. 
Of course, the Higgs sector may be sensitive to any refinements designed to solve the 
"/i problem". 



6 



The CLEO Collaboration H has quoted the 95% CL range 1 x 10"^ < B{b 57) < 
4.2 X 10""^ from their study of inclusive flavour changing neutral current decays of 5-mesons. 
Such low energy measurements provide constraints on any particular model framework where 
virtual effects of new particles can contribute to the decay. We have used the recent analysis 
described in Ref. I^Sj to compute this branching fraction within the GMLESB framework. 
The result of this computation is illustrated in the A vs. tan/3 plane in Fig. 3 for (a) yU < 0, 
and (b) fi > 0. Again, we have fixed M = 500 TeV. This constraint excludes a substantial 
region of the parameter space for negative /i due to constructive interference amongst the 
various SUSY and SM loop contributions. The CLEO experiment poses no constraint for 
yU > 0, since in this case the various SUSY loops contributions interfere destructively. The 
rate for 6 — > 57 could be sensitive to the modifications of the Higgs sector arising from new 
dynamics included to generate /i dynamically. For this reason, in Sec. Ill we include in our 
study parameter space points for which this 6 — > 57 constraint is not satisfied. 

Before turning to phenomenology, we note that the conditions for electroweak symmetry 



breaking determine the weak scale B in terms of tan (3 (while /i is fixed by the value of 
Mz). This value of B can then be evolved to Bq, its value at the messenger scale M. Since 
Bq is not generated at one loop, we expect that it should be small within our framework. 
Contours of Bq are shown in the A vs. tan /3 plane in Fig. 4 for (a) /i < 0, and (b) n > 0. We 
see that for positive values of /i, Bq is always very large. On the other hand, there is a region 
of parameter space with /i < where Bq is close to zero. If we take the model literally, we 
would conclude that tan P is fixed to be between 20 and 30, depending on the value of A and 
/i < 0. While this might be interesting in itself, this conclusion would probably be altered 
by the addition of new interactions that would be necessary to generate fi dynamically. For 
this reason we will remain agnostic about tan /3 and sgn /j, in the remainder of this paper. 



C. The decay of the Lightest NeutraUno 

We have seen that below the scale M, the GMLESB model looks just like the minimal 
supersymmetric model with (correlated) soft supersymmetry breaking terms, together with a 
very light gravitino as the LSP. The NLSP, which is usually the Zi, thus decays via Zi 7G 
and also via Zi —>■ ZG and Zi — > HiG if these decays are kinematically allowed. Expressions 
for these decay rates are given in Ambrosanio et. al. and will not be repeated here. 

The branching fractions for the various decays are shown versus A in Fig. 5 for (a) tan (3 = 
2,/i < 0, (b) tan/5 = 2,/i > 0, (c) tan/5 = 10,/i < and (d) tan/5 = 10,/i > 0. The 
messenger scale has been fixed at our canonical choice M = 500 TeV. We see that the 
photonic branching fraction dominates for the entire range of A even though m^^ is as 

heavy as ~ 180 GeV for A ~ 140 TeV. This is a reflection of the fact that Zi ~ B. Since the 
zino component of B is suppressed relative to the photino component by tan^^, we expect 
that B{Zi GZ) is suppressed by a factor sin^ 6w- The remaining suppression comes from 
the strong [3^ suppression of these decays. The decay rate to is negligible because the 
Higgsino component of Zi is tiny. 

The decay length L = I3-^_^'J^_^ct of the neutralino is illustrated in Fig. 6 for (a) M = 
500 TeV with /i < 0, (6) M = 5000 TeV with /i < 0, (c) M = 500 TeV with /i > and 
(d) M = 5000 TeV with /i > 0. The curves are for three different values of 75 = — ^ = 1.5, 2 
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and 4, (from bottom to top). In this plot, we have fixed tan (3 = 2. We expect our results 
are insensitive to this choice as long as Zi ~ 5. We see that for our canonical choice 
M = 500 TeV, the decay length varies from a fraction of a millimeter to a few centimeters. 
Thus the neutralinos will decay inside the detector and the displaced vertices from which a 
high energy photon shower emerges could provide additional confirmation of this scenario. 
For the larger value of M shown, the decay length could be as large as several meters if 
A/M is small, so that the neutralino would decay outside the detector. In this case, the 
topological signatures would be similar to those in the minimal model: the vestiges of the 
GMLESB model would show up only via the sparticle mass patterns. Intermediate values 
of M could cause the Zi to mainly decay outside the electromagnetic calorimeter, or within 
the muon chamber. Whether these decays can be readily identified and/or the photon 
energy measured is an important experimental issue. It is interesting to note that while the 
sparticle mass scale provides a handle on A, a measurement of the decay length of Zi 
would directly yield information about the messenger scale, particularly if the composition 
of the Zi could be determined from other experiments. 



III. SIGNALS FROM THE GMLESB MODEL AT THE TEVATRON 

The patterns of sparticle masses, and hence, the cross sections for various sparticle 
processes can be quite different from expectations in, for instance, the mSUGRA framework. 
In order to compute these cross sections as well as to generate SUSY events at the Tevatron 
within the GMLESB model framework, we have interfaced the output for the various weak 
scale parameters as obtained by RGE, starting from the GMLESB boundary conditions 



with ISAJET [17|. We begin by showing in Fig. 7 the cross sections for various SUSY 



processes as a function of A for the same cases {a)-{d) shown in Fig. 5. Again, we have 



fixed M = 500 TeV. We use the CTEQ2L structure functions p5| for our computations. We 
show the cross sections for the dominant W1Z2 and WiWi processes separately, and group 
together the processes of slepton and sneutrino pair production in the figure. The curve 
labelled "0th." refers to other chargino and neutralino processes while "Assoc." refers to 
the production of a gluino or squark in association with a chargino or a neutralino. We note 
the following. 

• Over the complete range of A where the cross sections are potentially observable, 
W1Z2, WiWi and slepton/sneutrino production processes dominate. The production 
of gluinos and squarks is always subdominant. This is a refiection of the fact that 
gluinos and squarks are rather heavy even for the smallest allowed value of A. 

• The strongly interacting sparticles get rapidly heavier as A increases, so that their 
cross sections drop off faster for a fixed collider energy. For the same reason, the pro- 
duction cross section for electroweak sparticles falls off the slowest, with the associated 
production cross sections in between. 

• The cross section for the production of "other" charginos and neutralinos is significant 
for smaller values of A shown, particularly if /x > 0. Presumably, this is because n is 
not overwhelmingly large and gaugino-Higgsino mixing tends to reduce the sparticle 
masses. 
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A. Event Simulation 



For each set (A, tan/3, M, sgn/i), of input GMLESB parameters, our RGE program 
yields a set of weak scale SUSY parameters. We use these as an input to ISAJET to 
generate SUSY events at the Tevatron. Thus, for any set of GMLESB parameters ISAJET 
generates all 2 — > 2 SUSY processes (those mediated by s-channel Higgs production must be 
run separately) with appropriate cross sections, and decays all sparticles as in the minimal 
SUSY model. The decay Z\ with a branching fraction of 100% (this is an excellent 

approximation as can be seen from Fig. 5) is added to the ISAJET decay table. Gravitino 



decays of sparticles other than Zi are ignored [26 



To model the experimental conditions at the Tevatron, we use the toy calorimeter simu- 
lation package ISAPLT. We simulate calorimetry covering —4 < 77 < 4 with a cell size given 
by fSj] X A0 = 0.1 X 0.0875, and take the hadronic (electromagnetic) calorimeter resolution 
to be Q.l /\fE (0.15/v^)- Jets are defined as hadronic clusters with Et > 15 GeV within a 
cone of AR = \/ Arj'^ + = 0.7 with \rij\ < 3.5. Muons and electrons with Et > 7 GeV 
and Irjil < 2.5 are considered to be isolated if the visible hadronic Et within a cone of 
AR = 0.3 about the lepton direction is smaller than 5 GeV. We identify photons within 
\ri^\ < 1 if Et > 12 GeV, and consider them to be isolated if the additional Et within a 
cone of AR = 0.3 about the photon is less than 4 GeV. Moreover, we assume that a photon 
within the acceptance is detected with an efficency of 80% (100%) if its energy is smaller 
(greater) than 25 GeV. In our analysis, we neglect multiple scattering effects as well as any 
detector-dependent effects such as lepton, photon or jet misidentification. Finally, in our 
simulation, we have not incorporated the finite decay length of the Zi but assumed that the 
Zi decays at the production vertex. This will introduce some small error in the direction 
of the photons from the decays Zi 'jG for our choice of M = 500 TeV. Although the 
displacement of the Zi decay vertices should be properly included in a complete simulation, 
we see from Fig. 6 that over most of the parameter range for which we have performed our 
simulation, the decay length is a fraction of a centimetre so that the results we show below 
should not be qualitatively altered. 



B. Classification of Events and Topological Cross Sections 

We classify GMLESB signals at the Tevatron primarily by the number of isolated pho- 
tons, and then further separate them by their lepton content. We also distinguish events 
which do or do not contain jets. 

In addition to a global cut I^j, > 30 GeV, we require that every event must satisfy at 
least one of the following lepton, photon or jet requirements which are motivated by the 
need for a trigger: 

• U with pT{i) > 20 GeV, or 2i with pr(^) > 10 GeV; 

• two isolated photons; 

• two jets with Et > 30 GeV and > 40 GeV. 

The results of our calculations of the cross sections for various event topologies after 
imposing the cuts and "trigger requirements" described above are shown in Fig. 8-Fig. 11 
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as a function of the parameter A. We have once again fixed M = 500 TeV and shown these 
cross sections for tan/? = 2 and 10, and both signs of fi as in Fig. 5. For each choice of tan/3 
and /i we show the cross sections for events with (a) no identified photons, (6) I7, and (c) 27. 
The sohd hnes show the cross section for events with jets, while the dotted hnes show the 
cross sections for "clean" events free of jet activity. We have performed event simulations for 
values of A between 40-140 TeV in steps of 20 TeV and denoted the cross sections n-lepton 
events by the symbol n in the figures. The following remarks about Fig. 8-Fig. 11 are worth 
noting: 

1. For all the four combinations of tan/3 and sgn/x shown in these figures, we see that 
cross sections for events with jets dominate the cross sections for clean events. This is 
because over most of the parameter space, W1Z2 and WiWi are the dominant sparticle 
production mechanisms, and at least Wi typically has a large branching fraction for 
hadronic decays. For the smaller values of A shown in the figure, Z2 Hr is the only 
two body decay channel that is kinematically allowed, so that the leptonic decays of 
Z2 dominate. As A becomes larger, the decays Z2 — >■ ZZi and Z2 — >■ HiZi become 
accessible and dominate the decay to right-handed sleptons, so that Z2 then mainly 
decays via its hadronic mode. This also accounts for why the dotted lines in Fig. 8- 
Fig. 11 exhibit a steeper fall-off than their solid counterparts. 

2. Most of the dilepton plus multijet events contain opposite sign dileptons in this case 
because contributions to the signal from gg and gq production are subdominant (see 
Fig. 7) since gluinos and squarks tend to be heavy. 

3. We see that for the jetty event sample, the cross section for I7 events is larger by a 
factor 1.5-4 than the cross sections in the or 2 photon event samples. This is, of 
course, sensitive to our assumptions about the 7 acceptance < 1) and detection 
efficiency and could be different for Run II of the Tevatron. 

4. For the case of the clean events shown in the figures, we see that cross sections where 
both photons are observed tend to be larger than those where the photons escape 
detection. This is presumably because it is easier for the photons to satisfy the isolation 
requirements than in the case of jetty events. 

5. A comparison of the four figures shows that the various cross sections vary rather 
weakly with tan j3 but show slightly more sensitivity to sgn fi (when tan j3 is not large). 

We have not made an attempt to compute SM backgrounds to the SUSY event sample 
from the GMLESB model. Background levels for the zero photon sample in case (a) shown in 
these figures have been previously estimated [^], although not with precisely the same cuts. 
We surmise that the presence of additional isolated photons, and possibly, also the presence 
of up to two significantly displaced vertices (without charged tracks emerging from them) 
would reduce the physics backgrounds to negligible levels. Of course, a careful computation 
that includes the effects of the non-zero decay length of the NLSP should ultimately be 
carried out to ensure there are no unforeseen surprises. 

Assuming that the signal is indeed rate-limited, we estimate the SUSY reach of the 
Tevatron for an integrated luminosity of {i) 100 pb~^, corresponding to the size of the 
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Run I data sample per experiment, (m) 2 fb~^, the integrated luminosity expected to be 
accumulated after about two years of Main Injector operation at design luminosity, and 
finally, {in) 25 fb~^, the integrated luminosity that might optimistically be accumulated at 



the proposed TeV33 upgrade ||18[ of the Tevatron. For our estimate of the Run I and Run II 
reach, we take the 5 signal event level as our criterion for observability in any one channel, 
while for TeV33, we take the observability level to be 10 events, and show the corresponding 
cross sections by the horizontal dashed lines in the figure. We see that when the data from 
Run I is analysed, the CDF and DO experiments will be probing A values of 50-60 TeV. With 
the Main Injector, experiments at the Tevatron should be able to explore up to A ~ 100 TeV. 
If TeV33 is able to accumulate a data sample of 25 fb^^, then the reach should extend out 
to about 135 TeV . For comparison with earlier studies of the Tevatron reach, these reach 
numbers correspond to rrig values of ~ 450, 800 and 1100 GeV. In contrast to the mSUGRA 
case p7| , p!8| , p8| , |29[| , we find here that the reach via the clean channels is smaller than via 
jetty channels. It should, of course, be remembered that we are not directly probing such 
massive gluinos at the Tevatron, but obtaining the signal via the chargino and neutralino 
channels. Also, we remind the reader that our estimate of the reach may be somewhat 
optimistic since we have assumed that backgrounds are completely negligible: there could 
be important detector-dependent backgrounds that may not be ignorable. On the other 
hand, the reach could be even larger than our estimate if we sum up the expected signal in 
the many different channels. 

Before closing this discussion, we remark upon the various attempts in the literature 
1^^141 to account for the e"'"e~77 event by the CDF collaboration |jl5| within the GMLESB 



framework. We see from Fig. 8-Fig. 11 that while it is indeed possible to have a cross section 
of about 10 fb for clean dilepton plus two photon events (corresponding to ~ 1 event in the 
Run I data), this event should have been accompanied by at least an order of magnitude (and 
possibly, as many as fifty) times as many events in other channels. For this reason, we feel 
that this interpretation of the CDF event is unlikely at least within this minimal framework. 

We should mention, however, that it is possible to reduce the ratio of slepton to elec- 
troweak gaugino masses by increasing the number of 5 -|- 5 fields in the messenger sector. 
This cannot be larger than four if gauge couplings are to remain perturbative up to the 
GUT scale, but it is possible to arrange for Wi — > iii^, ui and Z2 —>■ ii, vv to be the only 
two body decays of the charginos and neutralinos if A is not large. In this case, the hadronic 



signals from Wi and Z2 production would be greatly reduced ||30[. It would be of interest 
to simulate such a scenario to see whether it is possible for the dominant signal to be in the 
channel. But we stress that it is necessary to check all other signals that are likely 
to be present before attributing the CDF event to sparticle production within the GMLESB 
picture. 



IV. SUMMARY AND CONCLUDING REMARKS 

During the last year or two, we have witnessed the emergence of a phenomenologically 
viable alternative to the minimal SUGRA model for analyses of supersymmetry. This new 
scenario is similar to SUGRA in that SUSY is dynamically broken in a hidden sector of the 
theory which does not couple to the known particles and their super-partners via SM gauge 
interactions. Effects of SUSY breaking are communicated to the known particles via SM 
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gauge interactions with a messenger sector which also couples to the hidden sector. The 
difference is that while gravity plays the role of the messenger within the SUGRA framework;, 
new sets of messenger fields are invoked in these novel scenarios. Thus, while the messenger 
scale is necessarily 0{Mpianck) for SUGRA models, this scale may be as small as a few 
hundred TeV within the novel GMLESB scenario. 

Prom a phenomenological point of view, the SUGRA and GMLESB scenarios differ 
in two crucial respects. First, the boundary conditions for the RGE that determines the 
weak scale parameters of the theory are different: in the mSUGRA case, we have universal 
parameters at a scale ~ Mqut-, while in the GMLESB case, sparticle masses, which are 
radiatively generated at the messenger scale are proportional to the SM gauge couplings 
squared. As a result, sparticles with the same gauge quantum numbers have the same 
masses (unless they have large Yukawa type interactions) so that flavour changing neutral 
currents are automatically suppressed. Second, unlike the mSUGRA framework where the 
lightest neutralino is the LSP and the gravitino decouples from particle physics, the gravitino 
is superhght within the GMLESB scenario and decays into a gravitino and a photon (or, 
possibly also a Higgs or Z boson). For sparticles in the mass range accessible at the Tevatron, 
only the photon decays of Zi are significant, so that every SUSY event contains two isolated 
hard photons (though these may not both be detected). 

In Section II we have set up the parameter space for the simplest GMLESB scenario 
and examined the extent that this has directly been probed by experiments at high energy 
colliders or indirectly via the effects of virtual sparticles on the flavour changing decay 
6 — 57. We have also shown that a measurement of the decay length of the Zi yields 
information about the messenger scale. 

The main purpose of this paper was to study the cross sections for various event topologies 
that should be accessible at the Tevatron within the minimal GMLESB picture. Towards 
this end we have interfaced the weak scale SUSY parameters obtained from these boundary 
conditions with ISAJET to obtain these cross sections. We believe that our calculations 
are the flrst semi-realistic simulations performed within this framework. Our main results 
are exhibited in Fig. 8-Fig. 11 for experimental conditions suitable at the Tevatron. We see 
that, unlike in SUGRA where multijet plus ^j, events form the dominant event topology, 
multijet plus = 0, 1, 2 plus 1 photon events are the major component of the SUSY cross 
section. Similar events with two isolated photons or zero photons, which have only a slightly 
smaller cross section, may also be present at observable levels even in the Run 1 data sample 
(and certainly at the Main Injector upgrade) if any observation in the single photon channel 
is to be attributed to the minimal GMLESB realization of SUSY. It should be kept in 
mind that the relative sizes for the 0:1:2 7 cross sections are sensitive to our assumptions 
about the photon acceptance and detection efficiency. We also see that the cross section for 
clean multileptons plus photon event topologies is below our level of detectability during the 
current run: a handful of such events may, however, be present in the CDF and DO Run I 
data samples if A is not too large. At the Main Injector, up to several tens of clean 77 
and 7 plus multiple lepton events may be present. While the reach extends out to about 
A ~ 50 — 60 TeV for Run I experiments, the Main Injector (TeV33) should be able to probe 
up to A ~ 100 TeV (135 TeV), which corresponds to m-g ~ 800 GeV (1 TeV)! 

How stable are our conclusions to model variations? We note that there is a small region 
of parameter space (the cross-hatched region in Fig. 2) where fi is the NLSR For parameters 
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in this region, all SUSY events will contain at least 2-4 r leptons but no photons in the final 
state. We note that the phenomenology may be sensitive to assumptions about the messenger 
sector. For instance, if instead of assuming that it contains a single vector multiplet of SU (5), 
if instead we assume it contains four such multiplets (coupling roughly the same way), then 
slepton and squark masses reduce by about half relative to the electroweak gaugino and 
gluino masses. This could have a significant impact on sparticle decay patterns and the 
resulting phenomenology. 

In summary, we have examined the implications for experiments at the Tevatron in 
a new class of models where SUSY is broken at relatively low energy, and the effects of 
SUSY breaking communicated by gauge interactions. The production of sparticles at the 
Tevatron would then result in a variety of events with n-jets + m-leptons +A;-photons + 
I^rp. We have computed the cross sections for these event topologies under experimental 
conditions appropriate to the Tevatron, and mapped out its reach within the parameter 
space of the model for both the current run as well as the Main Injector (and TeV33) 
upgrades. Observation of these events would not only be spectacular in that it would signal 
the discovery of a fundamental new symmetry of Nature, but also in that it would imply 
the existence of a whole new family of particles not very far beyond the multi-TeV scale. In 
contrast to the case of the desert hypothesis, we would have a hope of directly probing this 
sector in the foreseeable future. In the mean time, we might be able to obtain |31] indirect 
information about the physics of this sector via the experimental determination of masses 
and other properties of sparticles. 
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NOTE ADDED: While this paper was in preparation, we received two papers ||32|J33 
where the spectroscopy and the phenomenology of the minimal GMLESB model as well 
as of models with extended messenger sectors is considered. While parts of these papers 
overlap with the present work, neither of these papers performed explicit event generation 
for the GMLESB model for the Tevatron collider. 
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FIGURES 



FIG. 1. Renormalization group trajectories for the soft SUSY breaking scalar masses and the 
gaugino masses Mj versus renormahzation scale Q from the messenger scale (M = 500 TeV) to the 
weak scale. In this example, we take A = 40 TeV, tan/3 = 2, /j, <0 and nit = 175 GeV. 

FIG. 2. Contours of (solid), ii (dashed) and Zi (short-dashed) and Wi (dotted) masses in 
the A vs. tan/3 plane of the GMLESB model with a single 5 + 5 representation in the messenger 
sector for a) jd < and 6) /x > 0. We also show contours for rUg (solid), ruq (dashed), At 
(short-dashed) and /x for c) < 0, and d) /j, > 0. We take rrit = 175 GeV and fix M = 500 TeV. 
The bricked regions arc excluded by theoretical constraints discussed in the text, while the shaded 
regions are excluded by experiment. The cross-hatched region is where fi is the NLSP. 

FIG. 3. Contours of the branching fraction (xlO^) for the decay b in the A vs. tan/3 

plane for M = 500 TeV, for a) fi < 0, and b) n > 0. The CLEG experiment has measured this 
branching fraction to be between (1 — 4.2) x 10~^ at 95% CL. 

FIG. 4. Contours of the value of the B parameter as obtained from the conditions of electroweak 
symmetry breaking but evolved to the messenger scale taken to be 500 TeV for a) /x < and & /x > 0. 

If there are no new interactions, the value of this parameter should be small so that the model 
taken literally picks out /x < and tan /3 =20-30. See the text for a further discussion of this point. 

FIG. 5. Branching fractions of the lightest neutralino Zi decays to the gravitino versus A for 

a) tan/3 = 2, /i < 0, b) tan/3 = 2, // > 0, c) tan/3 = 10, /x < and d) tan/3 = 10, /x > 0. The 
messenger scale is fixed to be 500 TeV. 

FIG. 6. The decay length in centimeters of the lightest neutralino Zi decays to the gravitino 
versus A for a) tan/3 = 2, /x < 0, M = 500 TeV and 6) tan/3 = 2, /x < 0, M = 5000 TeV. The 
three curves from bottom to top for 7-^ = 1.5, 2 and 4. Frames c) and d) are identical to a) and 

b) above except that /x > 0. 

FIG. 7. Total production cross sections for sparticles in the GMLESB scenario for the Tevatron 
collider operating at ^/s = 2 TeV. We show frames for a) tan/3 = 2, /x < 0, 6) tan/3 = 2, /i > 0, 

c) tan/3 = 10, /x < and d) tan/3 = 10, /x > 0. "0th." refers to other chargino and neutralino 
processes and "Asso." refers to the production of a chargino/neutralino in association with a 
gluino/squark. 
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FIG. 8. Topological cross sections from sparticle production and decay versus A in the GMLESB 
framework for the Tevatron collider operating at -y/s = 2 TeV with cuts and trigger conditions 
listed in the text. The messenger scale is fixed to be 500 TeV. We take tan/3 = 2 and < 0. We 
show frames for a) events containing no isolated photons, b) events containing a single isolated 
photon, and c) events containing two isolated photons. The solid curves correspond to events 
containing jets, while the dotted curves correspond to clean topologies (no jets). The curves are 
labelled according to the number of isolated Icptons present in the signal. The dashed horizontal 
lines correspond to the approximate reach of the Fermilab Tevatron with 0.1, 2 and 25 fb~^ of 
integrated luminosity with observability criteria listed in the text. 

FIG. 9. The same as Fig. 8 except that tan/3 = 2, /x > 0. 

FIG. 10. The same as Fig. 8 except that tan/3 = 10, /j, < 0. 

FIG. 11. The same as Fig. 8 except that tan/3 = 10, /Lt > 0. 
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